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. We analyse the properties of a very simple "balls-in-boxes" model which 

„^ , can exhibit a phase transition between a fluid and a condensed phase, 

\Q I similar to behaviour encountered in models of random geometries in one, 

two and four dimensions. This model can be viewed as a generalisation 
of the backgammon model introduced by Ritort as an example of glassy 
\& ' behaviour without disorder. 
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Introduction 

•i-H - 

Given the success of the dynamical triangulation approach to simulating two 
dimensional quantum gravity [1J], there has been a concerted effort to employ 
similar methods in understanding higher dimensional discretized gravity mod- 
els [§. In three and four dimensions one does not have the underpinning of 
matrix models to give confidence in the simulation results, but the general 
phase diagram of such models is now fairly clear. 

One finds a two-phase structure with an elongated/branched polymer phase 
and a crumpled/collapsed phase. Recent results suggest strongly that the tran- 
sition, originally thought to be continuous, is first order [EJ. In addition, the 
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vertex order distribution in the crumpled phase displays some peculiarities 
In D-dimensional models a singular D — 2 simplex appears in the crum- 
pled phase and mops up an extensive fraction of the volume. 

A similar phase structure is also encountered in one and two dimensional 
systems |6|-|8| but only recently was it hypothesized that the same mechanism 
might be responsible in all the cases @. The suggestion in || was that an 
effective polymer model described not only the elongated phase but also the 
essential features of the transition and the crumpled phase. In |§ it was pointed 
out that the polymer model admitted a description in terms of balls in boxes, 
with a box being a vertex and the number of balls in a box being the vertex 
order. The balls in boxes metaphor has already occurred in the analysis of 
glassy dynamics fIQ |, where the name "Backgammon models" was coined. 



In this paper we analyse the density driven transition of || using the methods 
of [10]. We first discuss the balls-in-boxes model in general terms, before going 
on to describe the transition itself. In closing we highlight the similarities with 



Bose-Einstein ||iT| , p!2]| condensation and the treatment of the spherical model 
phase transition in the classical paper of Berlin and Kac ]n|. The mapping 
between the balls-in-boxes model and the BP model of [[!]] is made completely 
explicit in an appendix. 



The balls— in boxes model 



Consider the system of N positive integer numbers: i — l,...,iV. We 
define the partition function: 



Z{N,M)= Yl P(Vi)---P(<lN)8 qi + ... + qNj M- 



n,—,qN 



The p(q) is some given function on positive integers taking on non-negative 
real values. If it is interpreted as a probability we have the additional require- 
ment that YlqP(q) — 1- We can view this system as an ensemble of M balls 
distributed into N boxes. Each box is weighted by the factor depending only 
on the number of balls it contains. In the glassy dynamics applications the 
choice in |TD| was p(q) = exp(/?5 9)0 ) an d there was no static phase transition. 



The natural choice of p(q) for models associated with random geometries is 
p(q) ~ as such factors appear as discretizations of terms coming from 
the metric integration measure. In what follows we will keep our discussion as 
general as possible, specializing to p(q) ~ g _/3 where appropriate. 

We are interested in the system with a fixed density, i.e. the average number 
of balls per box, p = M/N. Inserting the integral representation for the delta 



2 



function into (|l[) we obtain: 



Z(N,p) = -L J2 Pin) ■ ■ -PM JdXe-^-^-pN) = 

— I dAe lA ^ f Ep(?) e " iAg J = ^ / dAex P (^V + logF(e~ lA )) (2) 
where F(z) = E^i p(?) z q ■ 

The generating function F(z) contains all the information about the system. 
In this paper we will consider only the case when this function has a finite, 
non-zero radius of convergence z c . Notice that without the loss of generality 
we can put z c = 1, since due to the constraint q\ + . . . + q^ = M the change 
p(q) ~* p{(l) z c q introduces only a constant factor z~ M in front of the sum 
defining the partition function (|l|). From here on we keep z c — 1. We further 
consider two subcases : (a) the derivative of F(z) at the radius of convergence: 
F'(X) — J2qQP(o) is finite, and (b) F'(l) is infinite. 

The integral (^) can be calculated by the saddle point method. By deforming 
the contour in the lower half plane Im A < we find that the saddle point is 
the solution of the equation : 



p=-±\o g F{e-»)=G^) (3) 
dp 

and lies on the imaginary halfaxis Ao = — «/x for a certain range of p. To see this 
we notice that the function G(p) increases monotonically between the values 
Cr(oo) = 1 to G(0) = F'(l)/F(l) when the argument p goes from infinity to 
zero. As we will see G(0) corresponds to the critical density and we denote it 
by p c . Since the function G(p) is monotonic the equation (|]) can be solved for 
real positive p = p{p) with p between 1 and p c . The leading part of the free 
energy density f(p) calculated by the saddle point method is then 



f(p) = ^\ogZ(N,p)=pp + \ogF(e~n (4) 
and grows from zero to 



f(p) = ^\ogZ(N,p)=logF(l) (5) 

when the density grows from one to p c . This value of the free energy is in fact 
the upper limit for the free energy density and corresponds to the situation 
when all, or almost all, the probabilities p(qi) are decoupled : 
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/ \ N 

Z(N,M)= p{Qi)'--p(Qn)= EK?) =e iVlo g F ( 1 ). (6) 

91 j— ,QN V 9 / 

When p is larger than p c the free energy is blocked at the limiting value 
log-F(l) and does not change further with p. Thus, there are two phases: for 
p < p c the free energy in the thermodynamic limit is given by @ and (f|), 
while for p > p c the free energy is given by (|5[). To distinguish those two phases 
one can define the following expression 



%=M (7) 

which vanishes above p c and can serve as an order parameter for this sys- 
tem. Notice that the terms proportional to the derivative dp /dp arising from 
differentiating (§J) vanish due to the saddle point equation (|]). 

We will see later that for p > p c the finite size corrections to the formula (|5]) 
are of the form log p((p — p c )N) and come from the energy of a singular box 
containing an extensive (p — p c )iV number of balls. This is the phenomenon of 
condensation, and we call this phase the condensed phase. The other phase, 
for p < p c , we call the fluid phase. 

In case (b) the situation is different, since F'(l)/F(l) is infinite, and therefore 
the critical density p c is also infinite. This means that the system never enters 
the condensed phase and has only the fluid one. Similarly one can check that 
the model has no condensed phase when the radius of convergence of F(z) is 
infinite. 

The expression for the free energy is identical to those obtained in the branched 
polymer (BP) model J7j which has p = 2. We show explicitly in the Appendix 
that the BP model can be mapped directly onto the balls-in-boxes model. 



Critical exponents 



In this section we will calculate critical exponents for the specific case : p(q) ~ 
q~P . As we have already noted this is the natural weight to consider for models 
associated with random geometries (and is that used in the BP model of [^]). 
These weights are also known to appear in the broad class of models with self 
organized criticality (see for example ]L4|] and references within). Adding an 
exponential prefactor to the weights p(q) ~ e Kq q~^ will not affect the phase 
structure of the model, since as discussed before, it results only in the factor 
e -nM j n f ron t f the partition function Z(N, M). 
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The generating function is thus given by the Dirichlet series [|15[ F(z) 
J2^Li q~ l3 z q , and the critical density is given by: 



p c = F'(l)/F(l) = MLJl (8) 

where ({(3) is the Riemann Zeta function. For (3 > 2 the critical density is 
finite. It diverges when (3 approaches two and the system never enters the 
condensed phase. For (3 smaller than two the model has only the fluid phase. 
When the density approaches the critical one Ap = p c — p — > + , then the 
free energy density becomes singular. The first derivative of the free energy 
density with respect to p is equal p(p) which is singular. To see this let us 
expand the function G(p) on the right hand side of (|3]) for small p. There are 
two cases. For 2 < (3 < 3 the expansion starts from the singular term M : 



Ap = G{p) = b^- 2 + ... (9) 

where b is a positive constant depending on f3. For (3 > 3 the dominating part 
comes from the analytic part of G{p) : 



Ap = G{p) = ap + ... + bp /3 ~ 2 + ... (10) 

where a is also a constant depending on (3. The first dots denote other powers 
of p n coming from the analytic part for n < (3 — 2. Inverting the last two 
equations to solve for /i gives, in the limit Ap — > + : 



Of Ap 1 /^" 2 ) for2</?<3 

T" = Z 1 ~ \ ( n ) 
d P { Ap^ for (3 > 3 

For integer (3 one has logarithmic corrections, in particular for (3 = 3, p = 
Ap log Ap. By similar reasoning one can find the singularities of the free energy 
with respect to the parameter (3 @. 



Dressed box-occupation distribution 



If it were not for the constraint the probability of finding q balls in one box 
would be equal to the bare probability p(q). Because of the constraint the rela- 
tion between the effective or dressed probability 7r(q) and the bare probability 
p(q) is given by: 
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n{q) = Z(N,M) ? ' ' - p ( qN)6 q + q 2 + --- + qN,M 



Q2—,qN 



. s Z(N-l,M-q) 
= ™ Z(N,M) <12) 

The above expression is simply derived by observing that when q is fixed in 
one box the other N — l boxes must contain M — q balls. In the thermodynamic 
limit p = const and N — > oo we can use the saddle point result 



Z{N, M) ~ exp (jjiM + log F(e^)iV) (13) 
to obtain: 



in the fluid phase and 



in the condensed phase. 

We can now repeat the same reasoning for the two (or multi)-point probabil- 
ities. We fix the occupation numbers q a and q b in two boxes and find that the 
dressed two-point probability is : 



( \ ( W , Z(N-2,M-q a -q b ) 
K{q a ,qb) =P{qa)p{q b ) ^ (16) 

Calculating the last expression in the thermodynamic limit and comparing 
it with the one point probabilities and finds that in this limit occupation 
numbers in different boxes become independent, namely : 



n(Qa, qb) - ^{q a )^{qb) -> o (17) 

and the same holds for other the higher multi-point probabilities. 
For every finite N the distribution ir(q) must satisfy the constraint: 

X>(9)=p. (18) 
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It is easy to check that this constraint is satisfied by the the limiting distri- 
bution fll4|) in the fluid phase and at the transition. Exactly at the transition 
where /i = the constraint is satisfied by the bare distribution p(q). As we 
move into the phase condensed phase p > p c , the constraint is not satisfied by 



(|T5|), which implies that the limit is not uniform. 



In the fluid phase the the constraint (|18|) can be enforced by the exponential 
term in (p^|). In the condensed phase this would require e -M = z > 1 which is 
impossible because of singularity. The probability (|^) stays equal to the bare 
probability and the system compensates otherwise for the "missing" particles. 
This change in behaviour produces a phase transition, and is perhaps best 
understood by looking at the behaviour of n(q) for a finite number of boxes. 



Finite number of boxes 



The method described in the previous sections works only in the thermo- 
dynamical limit. As pointed out already the limiting distribution (|TJ|) does 
not satisfy the constraint and indicates some substantial finite size effects. To 
study those effects we used two methods: 

(i) Exact calculation of partition function. 

The partition function ([I]) satisfies the following recursion relation: 

M 

Z(N,M) =Y,Z(N -1,M - q)p(q), N>2 (19) 

q=l 

with initial conditions Z(l,q) = p(q). 



Equation (|19|) can be used to calculate the partition function and in 
principle any of its derivatives. The cost of calculating Z is proportional 
to NM 2 but the actual limitation on the size of the system which one 
can study using this method is given by the numerical precision as the 
value of Z(N, M) drops exponentially with N. We have found that the 
practical limit is about iV = 500. 
(ii) Monte-Carlo simulations. 

For the bigger systems we are left so far with the necessity of MC sim- 
ulations. The construction of the appropriate program is standard and 
we will not describe it here. The program allows to simulate the sys- 
tem of few thousands boxes with few hours of CPU time on an average 
workstation. 

Using both these methods one confirms the following picture of the system: In 
the fluid phase the distribution nicely converges to the limiting distribution 
(IT3) as expected. In the condensed phase the distribution follows the limiting 
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1000 



Fig. 1. The effective one-point distributions for different p and 500 boxes. The solid 
line is the bare distribution p(q) ~ q~ 4 . The finite-size corrected ir(q) distribution 
for p = 1.05 is given by the dot-dashed line, while the long dashed line is the 
thermodynamical limit. 

power-like distribution ( |T5D for smaller q but develops a peak in the tail. These 
effects are clearly seen on the figure [1] which illustrates the phase transition 
produced by change in balls density. 

One can include the finite size effects into the formula (|bl| ) by noting that \i 
depends on p and hence on q since p = (M — q)/N. The predicted distribution 
is plotted on the figure |IJ (dashed line) and the agreement with the exact 
result is excellent. The same reasoning can be applied to the formula (p!5|). 
This results in sharp decrease of distribution at large q, clearly indicating a 
different nature of finite size effects in this case. 

The peak in the tail appears as the system tries to compensate for the particles 
missing from (|T^). We can approximately predict the position of this peak by 
the following reasoning: The number of missing particles is 



N(P - E =N(p- = N(p - p c ) (20) 

so assuming that all the missing particle are in one box this is the position 
of the peak, which as moves linearly with increasing N. In figure [2] we plot 
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Discussion 



The density driven transition obviously has a strong similarity to Bose-Einstein 
condensation 11]], especially as presented in, for instance, [Q. Although this 
is usually formulated in a grand canonical ensemble with a chemical poten- 
tial the discussion may be paraphrased into language very similar to that 
used above. If one fixes the number of energy levels to iV and the number of 
particles to M = pN the partition function for Bose particles is 



E 



N 

J|exp( 



N 



-/%e 4 ) 6Q2 qi -M) 



(21) 



qi,...,q N i=l 



i=l 



where each box has its own energy level e,. If one slips into a continuum 
language the expression one arrives at for the particle density (in three di- 
mensions) is just 



P 
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where p is now the particle density per unit volume, s is a spin degeneracy 
factor and F(z) = Y%Li q~ 3 ^ 2 z q . At fixed particle density per unit volume 
z is pushed to one as T is reduced, at which point a sizeable proportion 
of the particles begin to condense in the ground state. Thus at the formal 
level the mechanism for the Bose-Einstein condensation and the transition 
discussed in this paper are identical, namely, a Dirichlet series is pushed to 
the limit of its radius of convergence. The similarities with the treatment of 
the spherical model phase transition in |L3| are also apparent. In that case 
too the (spherical) constraint is exponentiated and a saddle point "sticking" 
signals the phase transition. 

We would like to stress, however, that in spite of formal similarity there exists 
a substantial difference. In the Bose-Einstein condensation case each box has 
its own energy level and the energy is linear in number of balls. The system 
condenses into the lowest state. In the Backgammon case the boxes are equiv- 
alent and condensation is a kind of symmetry breaking process. Moreover, the 
energy of the box is not linear in number of balls which is absolutely crucial 
for existence of phase transition, so there is no direct interpretation in terms 
of occupied states. One way of looking at the model would be to imagine a set 
of iV "atoms" and M "electrons" . The electrons would be bound to the atoms 
and the energy of an atom with n electrons would be log(n), corresponding to 
the energy levels E{k) = 1/k. 

Variations of the model with different vertex weights can also be easily accom- 
modated into the saddle point approach. For instance, changing the weight of 
a vertex of order q from l/q@ to t q /q@ gives the following equation for the 
critical density 



UJ- 1) • q' 



1^ 



C(0) + (- 



(23) 



For univalent vertices we can see from the above that as t\ — > the asymptote 
for p c moves from one to two, so we reproduce the result of that the 
transition is suppressed for finite (3 when p = 2. 

In the case of infinite radius of convergence of F(z) the condensation does not 
take place. Clearly, this will happen when we restrict the number of balls that 
are allowed in one box (p(q) = 0, q > gmax) as the function F(z) is then a 
polynomial in z. 

The appearance of the ( function in the saddle point equation also encourages 
one to play number-theoretic games with the vertex weights. It is possible 
to regard the standard ( function as a bosonic partition function and define 
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"arithmetic" k-parafermions |17] with the use of a generalized Mobius function 



fi k (q) = l if q = \[^) r % 0<r,<A;-l 



(24) 



where the pi are the prime factors of q. The case k = 2 corresponds to arith- 
metic fermions and k = oo to bosons (the standard £ function). If the vertex 
weights are taken to be Pk^/q 13 the saddle point equation becomes 



CG9-1) C(*0-1) 



CO?) 



(25) 



The sort of fermionic partition functions generated here are thus not useful 
for suppressing a collapse transition. The vertex weights do not suppress large 
numbers of particles in a box, only those numbers with k identical prime 
factors are not allowed. The behaviour of function of (3 is thus not 

much changed, apart from the asymptote for large /3 being reduced to zero. 
These arithmetic parafermion models thus crumple at a lower p c for a given 
(3. 

The backgammon models are perhaps best thought of as models for the transi- 
tion in 4D gravity at the level of baby universes and their connecting "necks" 
II, identifying the polymer structure with the baby universe tree. As was 
noted in ||, the 4D simulations have a fluctuating density, as (almost) fixing 
N4 still allows iV and hence p to vary, but such a smearing of the density 
should not affect the saddle point equations and the resulting phase structure. 
The appearance of models with p(q) ~ q" 13 is not, as we have noted above, 
restricted to those arising in the context of discretized gravity simulations. 



In the model described in [14| the power exponent is generated dynamically 



and is calculated by assuming that the distribution function is properly nor- 



malised. This corresponds exactly to our constraint (18). This constraint is 



satisfied by a pure power law distribution only at the transition. 

The backgammon models of JlO| were originally formulated to investigate 



glassy behaviour, and their dynamics display many spin-glass like features. 
Although the energy landscape of the polymer inspired backgammon models 
here is a long way from that of the "golf-course" landscape of the original 
backgammon models it might still be a worthwhile exercise to investigate the 
dynamics of the polymer backgammon models in the light of the results of 
fLOf . The dynamics may prove as interesting as the static properties discussed 
in this paper. 
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Appendix A 



In this appendix we describe in detail the mapping between the branched 
polymer model of and the balls-in-boxes model. 

Let n q be the number of boxes with q balls in it. These numbers must satisfy 
the constraints: 



= £> ff = M (26) 

q g 

The number of configurations with given n q 's is given by (we assume 0! = 1) 



(27) 



no!ni! • • • n^! 
so the partition function (ffl) can be written as: 



, TV' 

Z(N, M) = £ \Po°pTpT ■ ■ -p(n) nN (28) 

where the sum runs over all sequences of positive integers satisfying ( pB"D . 

The number of rooted, planted, planar trees with V\ vertices of degree one, v 2 
vertices of degree two, and so on is given by (see for example 



1 JV! 

(29) 



N v\\ • • ■ vn^- 
These numbers satisfy 
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$> g = iV, £>, = 2iV-l 

1 <? 



(30) 



We see that the weight of a branch polymer depends on the vertex orders in 
the same way as a balls-in-boxes configuration depend on the box occupation 
numbers Therefore the BP partition function can be also written in the 
form with the factor (|2T| ) substituted by fl2"9|). If we now forbid the empty 
boxes by putting p(0) = and set M = 2N — 1 (p = 2) those two functions 
will be equivalent. 
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